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Abstract 

Nested structure, which is non-random, controls cooperation dynamics and biodi- 
versity in plant-animal mutualistic networks. This structural pattern has been ex- 
plained in a static (non-growth) network models. However, evolutionary processes 
might also influence the formation of such a structural pattern. We thereby propose 
an evolving network model for plant-animal interactions and show that non-random 
patterns such as nested structure and heterogeneous connectivity are both quali- 
tatively and quantitatively predicted through simple evolutionary processes. This 
finding implies that network models can be simplified by considering evolutionary 
processes, and also that another explanation exists for the emergence of non-random 
patterns and might provide more comprehensible insights into the formation of 
plant-animal mutualistic networks from the evolutionary perspective. 
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1 Introduction 



Elucidation of basic design principles behind mutualistic networks is a central 
topic in theoretical ecology for understanding cooperation dynamics and bio- 
diversity. In particular, nested structure is a well-known non-random structural 
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pattern observed in plant-animal mutualistic networks ( iBascompte et al.l . 120031 ) . 
In this structure, pollinators (animals) of a certain plant form a subset of those 
of another plant in a hierarchical fashion. Such non-random patterns often 
strongly control dynamics of ecological systems. For example, nested structure 



is bel ieved to minimize competition and increase biodiversity ( IBascompte et al 



20091 ). In addition to nestedness, the heterogeneous connectivity, which is de- 



fined as the occurence of diversified patterns of interaction amo ng species (both 



plants and animals), is fou nd in such mutualistic networks (jJordano et al 



20031 : ISaavedra et al.l . l2009[ ). Together with nestedness, heterogeneous connec- 



tivity might influence the dynamics of ecological systems. 



To elucidate the origin of n on-random patt e rns su ch as nested structure and 
heterogeneous connectivity, ISaavedra et al.l ( 120091 ) recently proposed the bi- 
partite cooperation (BC) model inspire d by food-web models based on species 
traits and external factors (reviewed in Stouffer et al. ( 20051 )). Although this 
model agrees well with real plant-animal mutualistic networks, it is a static 
(non-growth) network model in which the number of species (plants and pol- 
linators) is fixed. The structure of model-generated networks is determined by 
three types of parameters drawn from exponential or beta distributions: for- 
aging traits (e.g., efficiency and morphology), reward traits (e.g., quantity and 
quality) and external factors such as environmental context (e.g., geographic 
and temporal variation). 



The structure of interactions between plants and animals may be affected 
by not only species tr aits and external factors but also evolutionary history 
( iRezende et al.l . l2007al ) . Thus far, few efforts have been made to theoretically 
explai n the development of in teractions between plants and animals. For ex- 
ample, iRezende et al.l (l2007b( ) showed that plant-animal mutualistic networks 
are qualitatively more nested by using phylogenetic information in addition 
to species traits. Futhermore, they repoted that the phylogenetic signal was 
reported to be significant but weak in plant-pollinator network s. In addition, 
by using evolving network models, iGuimaraes Jr. et al.l ( l2007l ) explained the 
origin of frequency distributions of the number of interactions per species fol- 
lowing a power law with an exponential truncation (i.e., heterogeneous connec- 
tivity). However, this previous model was less predictive than the BC model, is 
a static (non-growth) network model (ISaavedra et al.l . l2009l ). Therefore, these 
results suggest that the evolutionary effects in plant-animal mutualistic net- 
works are not accurate. 



Nevertheless, evolutionary history may also determine the structure of plant- 
animal mutualistic networks because species are believed to mutually inter- 
act via traits or functions acquired through evolution. Ecological systems are 
very complex, and we have not yet understood the whole picture of ecologi- 
cal systems. To understand such complex systems, we need to consider them 
from several perspectives. Therefore, it is important to develop good evolving 
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network models that consider the effect of evolutionary processes on network 
formation. In the present study, we present an evolving network model of plant- 
animal interactions that takes into account simple evolutionary processes and 
use this model to investigate the effect of evolutionary history on the forma- 
tion of plant-animal mutualistic networks. Numerical simulations showed that 
this simple evolving network model better reproduced nestedness and hetero- 
geneous connectivity by using real data on plant-animal mutualistic networks. 
This result implies a not negligible contribution of evolutionary history to 
the structure of networks. Furthermore, we discuss the comprehensible ori- 
gin of nested structures and heterogeneous connectivity from an evolutionary 
perspective. 



2 The Model 



We developed an ev o lving network model based on the model proposed by 
Takemoto and Arital ( 2009 ) to explain plant evolution. In our model, a small 



initial plant-animal mutualistic network (Fig. [T] A) is preparated and repre- 
sented as a complete bipartite graph with hq plants and no animals, where 
uq is the initial number of species, and is an integer > 1. This initial network 
evolves according to two simple evolutionary mechanisms as follows: 

(i) Animals interacting with new plants are inherited from those interacting 
with existing (ancestral) plants. This mec hanism is based on the phyloge- 



netic resemblance between related species ( iGarland Jr. et al.l . 120051 ). We as- 
sume that new plants emerge because of mutations in ancestral plants. In 
our model, mutations in plants occur at the probability p and time t, and 
new plants originate from randomly selected existing plants. In this case, the 
traits (e.g., forms, colors, and chemical compounds) of new plants are similar 
to those of ancestral plants because of the phylogenetic resemblance. Since 
plant-animal interactions are based on these traits, the animals connecting to 
the new species are inherited from those connecting to ancestral plants (Fig. [1] 
B). By considering divergence, however, we model each animal such that it is 
inherited from animals interacting with the ancestral plant at the probability 
q (Fig. [1] C). This indicates that plants lose some traits due to the diver- 
gence. As explained above, since plants interact to animals via their traits, 
the plants also lose some interactions with animals. In addition, a bipartite 
network model gener ated based on the above inheritance (or copy) mecha- 



nism was proposed in lNacher et al.l (120091 ) to describe the evolution of protein 



domain networks at around the same time that we developed our model. 

(ii) Plants acquire new interacting animals after obtaining new traits (pheno- 
types) because plants are believed to connect to animals via phenotypic infor- 
mation as above. In particular, chemical compounds, which are often flower 
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pigment s, are involved in attract i ng ins ects (animals), protecting plants from 



insects ( iHarborne and Williams! . |2000| ). During evolution, living organisms 



(plants) are believed to obtain new phenotypes (traits) via genetic changes 
such as gene duplications and horizontal gene transfers. Since new phenotypes 
(traits) are based on existing phenotypes, plants interact with new animals. 
We modeled the above event to occur at the probability 1 — p and time t; 
a plant-animal pair is selected at random (Fig. [1] D) because new traits are 
based on existing traits, and thus, a plant acquires a new animal (Fig. [T]E). 

Our model has two parameters p and q. We can generate the model network by 
estimating these parameters using observable parameters in real plant-animal 
mutualistic networks: the number of plants P, the number of animals A, and 
the number of interactions L. 

Parameter p is estimated as 



where P = pt and A = [1 — p)t in our model. 

To estimate parameter g, we n eed to consider the time ev olution of L. This is 
derived as L {l—p)t/{l — q) (ITakemoto and Arital . 120091 ). Since A = (1— p)t, 
parameter q is estimated as 



1-^ 

L 



(2) 



Using Eqs. ([T]) and ([2]), we estimated parameters p and q from real data and 
generated corresponding model networks for comparison with real ones. 



3 Results and Discussion 



We first show that our model-generated networks reproduce nested struc- 
tures by using real d ata on plant-animal mutualistic networks reported by 



Saavedra et al.l (120091 ). The nested structure indicates that animals (pollina- 



tor) of a certain plant form a subset of those of another plant. The degree 
of nested structure is measured by nestedness A^, ranging from perfect non- 
nestedness (A^ = 0) to perfect nestedness (A^ = 1). When A^ is close to 1, we 
can see such a strong containment relationship in a hierarchical fashion. When 
A^ is close to 0, in contrast, plants have independent sets of animals: there is no 
such a relationship. In order to measure A^ of model and real-world networks. 
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Fig. 1. Schematic diagram of our model. The squares and circles represent plants 
and animals, respectively. (A) Initial plant-animal network is represented as a com- 
plete bipartite network. (B and C) The addition of a new plant and the resulting 
plant-animal interactions. The gray square represents a randomly selected plant. 
The filled square indicates a new plant resulting from duplication of the selected 
plant. The dashed lines represent possible interactions between the new plant and 
animals. (D and E) The addition of a new animal and the resultant interaction. The 
thick link between the gray nodes corresponds to a randomly selected plant-animal 
pair. The filled circle indicates a new animal. 



we u sed the BINMATNEST program (iRodriguez-Girones and Santamaria 
2006h . 



We also a ssessed the BC model an d random networks generated using the null 
model II (IBascompte et al.l . |2003| ) for comparison. The null model generates 
random bipartite networks with degree sequences that are similar to those of 
real networks. In the null model, the probability that a plant connects to an 
animal is proportional to the product of node degrees for a given plant and 
animal. 



The comparison of nestedness between models and real data is shown in Fig. O 
Our model and the BC model produced comparable nestedness in real network 
data. The null model could not predict nestedness. This result indicates that 
nested structure is independent from degree sequences and that the predictions 
of our model and the BC model are nontrivial. 



Furthermore, we closely evaluated the prediction accuracy of our model and 
the BC model by using the Pearson correlation coefficient and root mean 
square error between predicted nestedness Xj and observed nestedness i/i, de- 
fined as 
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Fig. 2. Comparison of nestedness between models and real data. The dashed line 
represents the perfect agreement between predicted nestedness and observed nest- 
edness. The values of nestedness from models are averaged over 100 realizations. 

where n is the number of samples. As shown in Table [H our model had higher 
prediction accuracy than the BC model. 

Table 1 

Comparison of prediction accuracy between our model and the BC model. 





Correlation Coefficient 


Root Mean Square Error 


Our model 


0.934 


0.0232 


BC model 


0.790 


0.0415 



Our model can clearly explain the origin of nested structure in plant-animal 
mutualistic networks by using fewer parameters. In this model, animals inter- 
acting with a new plant are inherited from those interacting with an ancestral 
plant because the two types of plants tend to be similar. New plants acquire 
some of the interacting animals via divergence (elimination of interactions). 
As a result, animals of an offspring plant become a subset of those of the 
parent plant, and thereby resulting in nestedness. 

Next, we investigated the frequency distribution of the number of interactions 
per species (degree distribution). In plant-animal mutualistic networks, we of- 
ten observe heterogeneous degree distributions, which follow a power law with 
an ex ponential truncation. Our na odel also reproduces heterogeneous distribu- 
tions ( iTakemoto and Arital . l2009h . which are not influenced by the structure 



of initial networks (i.e., no) in the case of large P and large A. The cumula- 
tive degree distributions of real plant-animal mutualistic networks (symbols) 
and model-generated networks (lines) are shown in Fig. [31 The cumulative 
distributions of only three representative networks are shown owing to space 
limitations. We could observe cumulative distributions of two types [Pcumikp) 
and Pcum(^A), where kp and /c^ denote the degrees on nodes corresponding 
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to plants and animals (pollinators), respectively] because plant-animal mutu- 
alistic networks are represented as bipartite graphs. The models studied were 
in good agreement with real data (Fig. 



(A) Marsh {Japan) (B) Grassland {Cass, New Zealand) (C) Subalpine (Japan) 




10° 10^ 10^ 10° 10^ 10^ 10° 10^ 10^ 

kft kft 



Fig. 3. Cumulative degree distributions in plant-animal mutualistic networks for 
plants (top column) and pollinators (animals) (bottom column). The circles and 
lines represent real data and models, respectively. The cumulative distributions of 
models are averaged over 100 realizations. 

In order to quantitatively verify goodness of fits between models and real 
data, we calculated the Kolmogorov-Smirnov (KS) statistics (distance) be- 
tween empirical distributions and predicted distributions for Pcum(fcp) (KSp) 
and Pcum(^A) (KS^). We also investigated the relationship between network 
size and the KS distance of degree distributions. Here, the network size and 
prediction accuracy were defined as P + A and the sum of the KS distances 
for the two distributions (i.e. KSp + KSa), respectively. 

The KS distances in our model and the BC model are compared in Fig. H] A. 
Overall, the KS distances were small for both models, i.e., in good agreement 
with real data. However, we also observed some low-accuracy predictions in 
both cases for small-sized networks. The KS distance increased as the network 
size decreased (Fig. |1]B) because the structure of small model-generated net- 
works may have strong fluctuated owing to its stochastic nature. The accuracy 
of our model was lower than that of the BC model for small networks. This 
might be because the BC model considers more parameters such as species 
traits and external factors than our model dose and is more independent from 
initial networks. In the case of large networks, however, our model more pre- 
cisely reproduced the degree distributions than did the BC model. 



Using our model, we could elucidate the origin of heterogenous distributions, 
which generally corresponds to mu ltiplicative processes or preferential attach- 
ments (IBarabasi and Oltvail . 120041 ). in plant-animal mutualistic networks from 
an evolutionary perspective. When new animals emerge, plants establish new 



7 



(A) 



□ KSp 



0.3 



E 
O 















• ..-ti 


' fl 




• 


• 


•• • 



^ • • • — 

0.1 0.2 0.3 

Our model 



(B) 



0.5 - 



I 



°-4 'r: a 



0.3 - oo 

o 

0.2 



0.7 I : : : , : : , : — 

Our model 

Q g _ BC model 



400 800 1200 1600 2000 

P+A 



Fig. 4. Kolmogorov-Smirnov (KS) statistics (distances) between models and real 
data. (A) Comparison of KS distance between our model and the BC model. The 
dashed line indicates that our model and the BC model had the same prediction 
accuracy. KSp and KSa are KS distances between observed and predicted distri- 
butions for -Pcum(^p) and Pcum(^A), respectively. (B) Correlation between network 
size and KS distance. 



interactions with them. Plants that interact many animals are likely to ac- 
quire more new animals because plant-animal pairs are selected at random. 
This phenomenon reflects preferential attac hment, which is similar t o that in 
the Dorogovtsev-Mendes-Samukhin model (jDorogovtsev et al.l . 120011 ). A sim- 
ilar mechanism exists in the case of animals. In our model, new plants inherit 
the animals of randomly selected plants. That is, animals also established new 
interaction with plants when new plants emerge because of random mutation. 
Thus, animals that interact with many plants tend to acquire more new plants. 
This is similar t o preferential attachment in the duplication-divergence model 
(jVazquezl . [20031 ). 



Our model is simpler than the BC model. In the BC model, parameters such 
as species traits and external factors have to be expressly provided; these pa- 
rameters need to drawn from given probabilistic distributions. Furthermore, 
the interactions between plants and animals are determined according to the 
linking mechanisms based on such parameters. In contrast, species traits and 
external factors and linking mechanisms are masked to a certain extent in our 
model, which takes into account evolutionary processes. Therefore, the above 
parameters, probabilistic distributions, and linking mechanisms based on pa- 
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rameters need not be directly considered in our model. However, this implies 
that species traits and external factors are not inessential. Species traits and 
external factors are important for plant-animal interactions. In our model, 
plant-animal interactions are determined as a result of species traits and ex- 
ternal factors. However, our model has no parameters corresponding to these 
traits and factors. Thus, it does not reveal the relationship of species traits 
and external factors with network structure in exchange for the simplification. 
This is a limitation of our model. 



We based our model on plant evolution. As explained in Sec. 2, plant chemical 
comp ounds are critical for plant-animal interactions ( iHarborne and Williamsl . 



2OOOI ) and hence, such interactions are often explained from the perspective 



of plants. Thus, we assumed that interactions between plants and animals are 
influenced by the traits of plants. Nonetheless, network formation may also 
be influenced by animal evolution. In this case, however, a more complicated 
model is required (i.e., more parameters). Therefore, we selected a model based 
on plant evolution. Due to the simplification, our model does not clearly ex- 
plain animals. For example, it does not consider what kind of animals join 
in networks when new species emerge. In other wards, it does not distinguish 
whether the new animals are related to existing ones or not, and it is not clear 
what traits animals have. For this reason, the strength of phylogenetic signal 
is not evaluated. Therefore, traits of animals and phylogenetic relationships in 
plant-animal mutualistic networks are not discussed by using our model. This 
is also a limitation of our model, and it is a problem needing to be solved in 
the future. 



It is surprising that the same model can explain two different systems: plant- 
animal mutu alistic networks and specie s-metabolite relationship in regards 
to fiavonoids (ITakemoto and Arital . l2009[ ). A hidden evolutionary relationship 
between fiavonoids and animals, through a radical speculation, may exist. 
Some fiavonoids, such as flower pigments, play roles in attracting insects (pol- 
linators) and protecting plants from insects. Such chemical compounds may 
correspond to traits (or phenotypes) of plants. Thus, the structure of species- 
metabolite relationships is expected be similar to that of plant-animal net- 
works. 



As shown above, structural patterns are also reprod uced through simple evolu- 
tionary processes. However, iRezende et al.l ( l2007bl ) showed that the phyloge- 
netic signal was significant but weak in an observed plant-animal mutualistic 
network. Our model, which shows high phylogenetic signal, does not explain 
this network well. This may be cause that there are several independent mecha- 
nisms for generating nested structure and heterogeneous connectivity. Network 
formation may be explained several independent theories. One is described by 
the BC model, which implies low phylogenetic signal. The other is based on 
our model, which suggests high phylogenetic signal. In nature system, the 
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nested structure and heterogeneous connectivity may be independently gen- 
erated by both mechanisms. Due to the balance of mechanisms, thus, we can 
expect to observe plant-animal mutualistic networks with high phylogenetic 
signal in addition to those with low phylogenetic signal. We believe that there 
are other plant-animal mutualistic networks descried by our model. Species 
traits, which are important for plant-animal interactions, are inherited from 
ancestral species to new species. Thus, evolutionary history may greatly affect 
the formation of plant-animal mutualistic networks. However, earlier models 
do not reveal the effect of evolutionary history on ecological networks because 
they do not clearly consider evolution. We believe that our model (an evolving 
network model) is important to obtain a deeper understanding of the forma- 
tion of plant-animal mutualistic networks. 

Regarding external factors, our model has a limitation. Heterogeneous envi- 
ronments (e.g., geographic, temperature, and nutrient differences) may reduce 
phylogenetic signals in plant-animal mutualistic networks. Our model may not 
clearly explain the formation of networks under such situations because it has 
no parameter corresponding to external factors due to the simplification. Our 
model can be applied to ecological networks if plant-animal relationships are 
comprehensively obtained under ideal conditions (e.g. environmentally homo- 
geneous islands). 



For simplicity, we did not consider a number of important evolutionary pro- 
cesses such as extinction and convergent evolution. In particular, degree distri- 
butio ns might be altered because of such extinctions (lEnemark and Sneppenl . 



20071 ). However, our results indicate that such mechanisms might have only 



negligible effects on the structural properties of ecological networks such as 
nested structure and heterogeneous connectivity. In plant species, genome dou- 
bling (polyploidity) is a major driving force for increasing genome size and 
the number of genes ( lAdams and Wendell . l2005l ). Duplicated genes typically 
diversify in function, and some are phenotypically expressed as forms, colors, 
syntheses of chemical compounds, etc. Such phenotypes (traits), especially the 
syntheses of chemical compounds, are therefore expected to increase, indicat- 
ing that the effect of node losses may be disregarded when considering the 
global tendencies of plant-animal mutualistic networks. However, this implies 
that deletions of nodes and interactions are not unnecessary. Such evolution- 
ary mechanisms might play important roles in determining partial (or local) 
interaction patterns of bipartite relationships. Thus, we need to focus on such 
evolution processes in the future to fully understand the formation of plant- 
animal mutualistic networks. 



Although our model has the above limitations due to simplification of model, 
it may be extensible by considering several factors (parameters). We believe 
that our model helps to reveal the relationships such factors and evolutionary 
history in plant-animal mutualistic networks and to more deeply understand 
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the formation of networks in the future. 



4 Conclusion 



By considering simple evolutionary processes, we have shown that non-random 
structural properties of plant-animal mutualistic networks, such as nested 
structures and heterogeneous connectivity, can be both qualitatively and quan- 
titatively predicted. This finding suggests that network models can be sim- 
plified by considering evolutionary processes. Furthermore, it implies that an- 
other explanation exists for the e mergence of no n -rando m structures (i.e.. 



from the evolutionary perspective). iRezende et al.l (j2007al Jb[) and our simple 



description of interactions between plants and animals indicate the impor- 
tance of evolutionary history, which is worth further investigation. The evolv- 
ing network model does not contradict static (non-growth) network models 
because these models describe respective different mechanisms for formation 
of non-random structures, and therefore provides additional insights into the 
formation of plant-animal mutualistic networks. 
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